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ENSO dynamics: low-dimensional-chaotic or stochastic? 

T. Zivkovic, 1 and K. Rypdal, 2 

Abstract. We apply a test for low-dimensional, deterministic dynamics to the Nino 3 
time series for the El Nino Southern Oscillation (ENSO). The test is negative, indicat- 
ing that the dynamics is high-dimensional/stochastic. However, application of stochas- 
tic forcing to a time-delay equation for equatorial-wave dynamics can reproduce this stochas- 
tic dynamics and other important aspects of ENSO. Without such stochastic forcing this 
model yields low-dimensional, deterministic dynamics, hence these results emphasize the 
importance of the stochastic nature of the atmosphere-ocean interaction in low-dimensional 
models of ENSO. 



1. Introduction 

Prediction is the ultimate goal of meteorology and cli- 
matology and the issue of predictability is crucial. Predic- 
tion in these disciplines are mostly probabilistic, but there 
may be different rationales for a probabilistic description. 
These rationales are intimately linked to the various mean- 
ings of the concepts of determinism, deterministic chaos, 
and stochasticity. The fundamental laws of classical physics 
are deterministic, since the future and past of the state of 
such a system is uniquely given by the state at a given time. 
The majority of models of weather and climate are deter- 
ministic in this sense. After the discovery of deterministic 
chaos, it has been realized that the deterministic evolution 
may be sensitive to the initial conditions, in the sense that 
small perturbations in the initial conditions grow exponen- 
tially with time. Predictability in such systems is limited 
by this exponential growth rate, given by the largest pos- 
itive Lyapunov exponent of the system. From a practical 
viewpoint the most important aspect of this insight is that 
predictability can be limited in this way, and may require 
a probabilistic description, not only in systems with a large 
number of degrees of freedom, but also in in simple, low- 
dimensional, nonlinear systems. 

Before the advent of chaos theory, unpredictability was 
considered to be a practical consequence of our inability 
to specify and solve the evolution equations for the micro- 
scopic state of high-dimensional systems. However, these 
days high-speed computing allows us, not only to solve 
numerically low-dimensional, nonlinear problems, but also 
high-dimensional general circulation models of the climate 
system. Nevertheless, the problem of sensitivity to ini- 
tial condition seems to persist, at least when it comes 
to weather prediction and climate variability on interan- 
nual to decadal time scales. The fluctuations of macro- 
scopic variables around some time-averaged or ensemble- 
averaged state could be described in a probabilistic man- 
ner by representing them as stochastic processes or alter- 
natively as low-dimensional chaotic processes. Thus, when 
dealing with data in climatology, either from observation 
or from large-scale simulations, the climate dynamicist will 
have to ask the question of whether prediction for the 
system/phenomenon of interest is better served by a low- 
dimensional chaotic model of climate variability, or by a 
high-dimensional (stochastic) model. 
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The answer to this question depends on whether the sys- 
tem/phenomenon represents a self-organization of the dy- 
namics into an effective small number of degrees of freedom. 
One approach is of course to reduce the high-dimensional 
model to a low-dimensional one via a series of approxima- 
tions and simplifications, i.e., via theoretical model reduc- 
tion. In many cases, this is very demanding and many differ- 
ent reductions are possible, and it is hard to know whether 
a reduced model is what Nature abides to. Hence, there is a 
demand for methods by which it is possible to decide from 
the observation data whether the system dynamics can be 
uniquely projected onto an attractor in a low-dimensional 
phase space. Such methods exist, based on Taken's time- 
delay embedding theorem [Takens, 1981], and from 1980 
onwards there is a growing literature on reconstruction of 
chaotic attractors from time series and computation of at- 
tractor dimension and largest Lyapunov exponent [Abar- 
banel, 1996]. The idea is to assume that the system state 
vector evolves according to the system of first-order ordinary 
differential equations describing the trajectory z(t) on the d- 
dimensional attractor in ap > d-dimensional phase space. If 
the system is autonomous and the attractor of the trajectory 
has dimension d, Takens' time-delay method [Takens, 1981] 
can be used to construct an m > 2d-dimensional embedding 
space on which the attractor can be mapped continuously 
and one-to-one. In practice this method works only if the 
attractor dimension d is reasonably low. Dynamical systems 
with a large number of independent or weakly dependent de- 
grees of freedom can only be described either by large-scale 
numerical simulation or by stochastic methods. For such 
systems the phase-space attractor is also high-dimensional 
and cannot be mapped one-to-one onto a low-dimensional 
time-delay embedding space. The computation of attractor 
dimension then typically fails to converge when embedding 
dimension m is increased, but such convergence can be dif- 
ficult to detect if the time series is short. Kaplan and Glass 
[1992, 1993] devised a direct test for the existence of low- 
dimensional deterministic dynamics which is useful for short 
time series. This is the kind of test that will be employed 
in this paper to the instrumental time series for the El-Nino 
Southern Oscillation (ENSO). 

The remainder of the paper is structured as follows: In 
section (2) we briefly describe a time-delay equation for 
equatorial wave dynamics and ENSO, and in section (3) 
and (4) we review the time-delay phase space reconstruction 
technique and the test for determinism. In section (5) these 
techniques are employed to the Nino data and to numerical 
solutions to the time-delay equation with and without sea- 
sonality subtracted and with and without stochastic forcing. 
We also apply a superposed-epoch analysis to these data to 
highlight the characteristic waveforms of ENSO episodes as 
manifested in the Nino signal. The implications of our find- 
ings are also discussed in this section and summarized in 
section (6). 



X - 2 



TATJANA ZIVKOVIC, KRISTOFFER RYPDAL: EL-NINO DYNAMICS 



2. A time-delay equation for ENSO 

An El-Nino episode is characterized by an increase in the 
sea-surface temperature (SST) in the eastern Pacific, with 
a strong impact at the coast of Peru. Generally, the SST 
in this area increases during the winter period, but occa- 
sionally (often every three to seven years) the temperature 
increase is more pronounced and this phase is called El- 
Niho. The opposite phase, a strong decrease in the SST, is 
coined La Nina. The intensity of El-Nino is higher than the 
intensity of La-Nina, and the SST distribution is positively 
skewed [Kondrashov et al, 2005]. The increase in SST is 
followed by stronger precipitation, which is distributed over 
the entire Pacific basin, while during La Nina and normal 
conditions, the precipitation is heavier in the western Pa- 
cific. This is also due to the Southern Oscillation (SO) in 
the atmosphere, which is characterized by lower pressure 
over the western Pacific during La Nina and normal condi- 
tions. According to the Bjerknes hypothesis [Neelin, 2011], 
there is an initial warming in the SST in the eastern Pa- 
cific, which weakens the trade winds from the east and gives 
rise to westerly wind anomalies, which further increase SST. 
This positive feedback (which is the El- Nino phase of ENSO) 
is also accompanied by the formation of the eastward Kelvin 
wave. Concurrently, a westward Rossby wave is formed in 
the middle of the ocean basin, and after it is reflected at the 
western boundary of the ocean basin, an eastward Kelvin 
wave is formed, which has a cooling effect. This wave causes 
the La Nina phase of ENSO. 

A large body of literature is concerned with the dynam- 
ical modeling of ENSO (a review can be found in Dijk- 
stra [2010]). An interesting class of models is described 
as delay-differential equations, which are linear and au- 
tonomous [Battisti and Hirst, 1989], nonlinear, autonomous 
[Suarez and Schopf, 1988], or nonlinear, periodically forced 
equations [Tziperman et al, 1994; Miinnich et al., 1990]. 
The latter will be our focus in this paper and can be formu- 
lated as an time-delay equation for the thermocline depth 
h(t): 



and that an observed time series x(t) is generated by the 
measurement function g : TV — > 7Z: 
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where k is an ocean-atmosphere coupling parameter, ck is 
the velocity of the wind-forced Kelvin mode, L is the ocean 
basin width, and cr is the velocity of the Rossby wave. Q. 
is the frequency of the seasonal cycle, while a, d and c are 
constants. The cosine function in equation (1) accounts for 
the annual periodicity in the SST data. This delay differen- 
tial equation has two time delays: n = L/ck + L/2cr, and 
T2 = L/2ck- Here, ri is the summation over a time it takes 
a Rossby wave to travel from the middle of the ocean basin 
to the western boundary and then be reflected as a Kelvin 
wave, while T2 is the transit time for the Kelvin wave which 
travels from the middle of the basin to reach the eastern 
Pacific. 

3. Phase space reconstruction 

Before we describe the test for determinism (or more 
precisely, for low-dimensional deterministic dynamics), we 
briefly explain how the phase space can be reconstructed 
from scalar time series of length N by time-delay embed- 
ding. Suppose that the phenomenon under study can be 
described by a state vector z(t) in a phase space of dimen- 
sion p, i.e., z evolves according to an autonomous system of 
1st order ordinary differential equations: 



c/z 

dt 



s(t) = ff («(t)). 



(3) 



Further, assume that the dynamics takes place on an invari- 
ant set (an attractor) A C 7Z P in phase space, and that this 
set has box-counting fractal dimension d. Since the dynami- 
cal system uniquely defines the entire phase-space trajectory 
once the state z(t) at a particular time t is given, we can 
define uniquely an m-dimensional measurement function, 



g : A -> TZ m , g(z) = (x(t), x(t + r), 



s(t + (m-l)T)), 



(4) 



where the vector components are given by equation (3), and 
t is a time delay of our choice. If the invariant set A is 
compact (closed and bounded), g is a smooth function and 
m > 2d, the map given by equation (4) is a topological em- 
bedding (a one-to-one continuous map) between A and lZ m . 
The condition m > 2d can be thought of as a condition for 
the image g(A) not to intersect itself, i.e. to avoid that two 
different states on the attractor A are mapped to the same 
point in the m-dimensional embedding space TV 71 . If such 
an embedding is achieved, the trajectory x(t) = g(z) (where 
g(z) is given by equation (4)) in the embedding space is a 
complete mathematical representation of the dynamics on 
the attractor. Note that the dimension p of the original 
phase space is irrelevant for the reconstruction of the em- 
bedding space. The important thing is the dimension d of 
the invariant set A on which the dynamics unfolds. 

4. Test for determinism 

This method was recently successfully applied to the stud- 
ies of magnetospheric organization during magnetospheric 
storms and substorms [Zivkovic and Rypdal, 2011, 2012]. 
When a system is low-dimensional deterministic, the direc- 
tion of the trajectory (its tangent) is a function of the po- 
sition in the reconstructed phase space (from equation 2). 
This means that trajectories emanating from points in a 
small neighborhood in phase space have almost parallel di- 
rections. On the other hand, corresponding trajectories in 
a stochastic or high-dimensional system have directions in 
a low-dimensional embedding space which are not uniquely 
dependent on the position in this space, and therefore the 
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Figure 1. A(r) averaged over results computed from an 
ensemble of ten realization of the O-U process (squares) , 
and for the realizations with randomized phases (trian- 
gles). Error bars denote standard errors of the mean A(r) 
(which is (l/\/l0)xthe standard deviation of the distri- 
bution of the ten samples). 
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Figure 2. The Nino 3 monthly time series. 



tangent can have a different direction the next time it re- 
curs to the same neighborhood. Let b denote a small time 
increment and envisage a portion of phase space spanning 
the entire attractor divided into an enumerable set of small 
"boxes" of size corresponding to the length of the trajectory 
increment: 



m-dimensional phase space is 



r _ n _ 1 ,2 1/2 F[(m + l)/2] 
n - Hn -^ [ ™' r(m/2) ; 



(8) 



where F is the gamma function. The deviation in (Vj) be- 
tween a given time series and a random-walk can be char- 
acterized by a single number given by the weighted average 
over all boxes of the quantity, 



A(r) 



(9) 



where we have explicitly highlighted that the averaged dis- 
placement (Vj)(r) of the trajectory in the reconstructed 
phase space depends on the time delay r. For a completely 
deterministic signal we have A(r) = 1, and for a completely 
random signal A(r) = 0, hence this quantity can be consid- 
ered as a measure of determinism. 

In Figure 1, we show A(r) averaged over ten numerical 
realizations of the Ornstein-Uhlenbeck (O-U) stochastic pro- 
cess, for embedding dimension is m — 8, and 6=1. The 



Ax(t) = [x(t + b) - x(t),x(t + r + b) - x(t + r), . . . , 

x(t + (m- l)r + 6) - x(t + (m- 1)t)}, (5) 

The tangent for the fc'th pass of the trajectory through box j 
is the unit vector Ufej = Axfc i3 (t)/|Axfc^(t)|. The estimated 
averaged displacement vector in the box is 




where rij is the number of passes of the trajectory through 
box j. If the embedding dimension is sufficiently high and 
in the limit of vanishingly small box size, the trajectory di- 
rections should be aligned and the length Vj = |Vj| = 1. In 
the case of deterministic dynamics and finite box size, Vj 
will not depend very much on the number of passes rij , and 
Vj will converge to 1 as rij — > oo. In contrast, for the tra- 
jectory of a random process, where the direction of the next 
step is completely independent of the past, Vj will decrease 
with rij as Vj ~ n a ■ The degree of determinism of the 
dynamics can be assessed by exploring the dependence of Vj 
on rij . In practice, this can be done by computation of the 
averaged displacement vector length: 

L n = {Vj) nj=n , (7) 

where the average is done over all boxes with same number 
n of trajectory passes. 

Let us recall that we are describing a test to distinguish 
signals described by low-dimensional dynamical systems, 
i.e., signals that in the continuous-time limit are solutions 
to differential equations and therefore continuous, and sig- 
nals described as a stochastic process. The archetype of 
a random and continuous stochastic process is the Wiener 
process (Brownian motion), which in discrete time is a ran- 
dom walk. The random walk has random increments, and 
hence the displacement vectors Ax(t) in an m-dimensional 
embedding space will have random directions. Hence, for a 
random walk L n ^oo — 0. For finite n, however, there will 
be a finite statistical spread of Vj, and as shown in Kaplan 
and Glass [1993], the average displacement of n passes in 
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Figure 3. (a): A(r) for the Nino 3 time series (squares) 
and A(r) averaged over ten different phase-randomized 
versions of this time series (triangles). The error bars 
of the latter denotes the standard deviation of A(r) over 
the distribution of these ten samples, (b): autocorrela- 
tion function r(r) for the Nino 3 time series. 
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Figure 4. Climatology for Nino 3 time series (red) and 
the sinusoidal climatology used in equation (1) (blue). 
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O-U process is described by the stochastic equation: 



dS t 



-XSt + adWt 



(10) 



where W(t) is the Wiener process. It is a more physically re- 
alistic random process for many phenomena than the Wiener 
process, since the damping term — XSt makes it bounded. 
In Figure 1 we also show mean A(r) computed from the 
same ten realizations, but after randomization of the phases 
of the Fourier coefficients. This randomization leaves the 
power spectral density, and hence the autocorrelation func- 
tion, unchanged. Hence, A(r) should also be unchanged for 
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Figure 5. (a): A(r) for original Nino 3 time-series with 
climatology subtracted (blue squares), and A(r) averaged 
over ten realizations of phase-randomized versions of this 
time series (red triangles). The error bars of the latter de- 
notes the standard deviation of A(r) over the distribution 
of these ten samples, (b): A(r) averaged over ten real- 
izations of time series generated by equation (1) with a 
white-noise stochastic forcing and with climatology sub- 
tracted (blue), and the same for ten phase-randomized 
versions (red). In (b) the error bars denote the standard 
deviation of the mean A(r). 
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a random process, which Figure 1 demonstrates. On the 
other hand, for a signal from a low-dimensional chaotic sys- 
tem, which has to be nonlinear to be chaotic, the random- 
ization of phases will destroy the nonlinear coupling between 
Fourier modes and make A(r) more similar to a random sig- 
nal, i.e., it will be reduced compared to A(r) for the original 
signal. Examples of this were shown in Kaplan and Glass 
[1992, 1993]. 

5. Results 

We analyze Nino 3 data (see Figure 2) obtained from the 
sea ice and SST data set (HadlSSTl) [Rayner et at, 2003]. 
Nino 3 data is the area-averaged SST from 5S-5N and 150W- 
90W, with monthly resolution and the time span between 
1871 and 2008. We compute A as a function of time-delay, 
which is also used in the phase-space reconstruction proce- 
dure. The embedding dimension is m = 10, 6 = 1, i.e., the 
box size is equal to the average distance between two succes- 
sive points on the reconstructed phase-space trajectory. It 
has been shown in Kaplan and Glass [1993] that increasing 
embedding dimension can increase A in deterministic sys- 
tems, while it should not influence A in stochastic systems. 
In principle, higher m is better, but we also have to con- 
sider the number of data points available for the test, and 
this decreases with increasing m (from equation 4). In our 
analysis, m = 10 seems a suitable choice of embedding di- 
mension. In Figure 3 we observe that A(r) exhibits spiky 
dips at those r where the autocorrelation function r(r) has 




Figure 7. A(r) averaged over results computed from 
ten simulated time series like that shown in Figure 6 (a) 
with climatology included (blue squares). After phase- 
randomization (red triangles). The error bars denote the 
standard deviation of the mean A(t). 




Figure 6. Time series analyzed in Figure 5.: (a): simu- 
lation, (b): Nino 3. 



Figure 8. Simulation of time-delay equation with no 
stochastic forcing, (a): n = 0.4088. (b): r a = 0.4198. 
(c): the signal in (b) de-seasonalized by wavelet filtering. 
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an extremum. This is a spurious feature of the technique 
which is explained in Kaplan and Glass [1993]. We also plot 
A(t) averaged over an ensemble of ten surrogate time series 
with randomized phases. 

From Figure 3, we can conclude that the dynamics 
underlying the Nino 3 data is dominated by a nonlin- 
ear, low-dimensional component, since A(r) for the phase- 
randomized time series is considerably reduced compared to 
that computed for the original Nino 3 data. However, as we 
will demonstrate next, the nonlinear and low-dimensional 
behavior is rather trivial and derives from the seasonal cy- 
cle. This cycle can be represented by the climatology, which 
is the mean over all data for every month of the year. In 
Figure 4 we plot the Nino 3 climatology along with the sinu- 
soidal climatology used in equation (1). Next, we subtract 
the climatology from the Nino 3 data, and compute A(r) 
again (Figure 5 (a)). A(r) is strongly reduced and random- 
ization of phases makes no discernible difference, indicating 
the the positive test for determinism in the original Nino 3 
time series is an effect of the seasonal cycle. 

This implies that the Nino 3 time series with the seasonal 
cycle subtracted can be described a stochastic process sim- 
ilar to the O-U process for which results were presented in 
Figure 1. It does not imply, however, that the time-delay 
equation (1) has to be abandoned, since inclusion of stochas- 
tic forcing term may give rise to such stochastic dynamics 
and produce results compatible with the test for determin- 
ism in Nino 3 data. In the following we shall first show that 
such a time-delay equation with stochastic forcing is able to 
produce time series with these properties. Next, we shall 
demonstrate that in absence of stochastic forcing, the time- 
delay model tends to produce low-dimensional, deterministic 
dynamics, even when seasonality is removed. 

In order to investigate the effect of such stochastic forc- 
ing we add a term aw(t), where w(t) is a Gaussian white 
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Figure 9. (a): A(r) averaged over results computed 
from ten simulated time series like that shown in Fig- 
ure 8(b) (blue squares). After phase-randomization (red 
triangles), (b): A(t) averaged over results computed 
from ten simulated time series like that shown in Fig- 
ure 8(c) (blue squares). After phase-randomization (red 
triangles). Error bars denote standard deviation of the 
ensemble mean. 



noise with unit variance, to equation (1) and solve the equa- 
tion numerically with the parameters a = 0, d = l,k = 
100, c = 1.3, a = 0.7 and ti = 0.2. Figure 5(b) shows 
A(t) ensemble-averaged over ten realizations of the simu- 
lation after the climatology has been subtracted. In the 
same plot we show A(t) averaged over the ten correspond- 
ing phase-randomized realizations. The low values of A(r), 
and their non-responsiveness to phase randomization, indi- 
cate that for these parameters the de-seasonalized solutions 
to the stochastically forced equation (1) have the character 
of a stochastic process similar to the de-seasonalized Nino 
3 time series. A qualitative similarity is also apparent from 
Figure 6, where we plot a realization of this process along 
with the Nino 3 time series, both with the climatology sub- 
tracted. Both time series have been normalized to have 
unit variance. This similarity is also apparent when A(r) 
is computed from the simulation without subtracting the 
seasonal cycle, as shown by comparing Figure 7 with Figure 
3. The model described above, but without stochastic forc- 
ing (a = 0), was studied by Ghil et al. [2008]. They used the 
same choice of parameters a, c and d as above, but varied 
the parameters k and n. They computed and plotted the 
global maximum of the thermocline depth h(t) as a function 
of ri and k, and observed discontinuities in this parameter 
space separating regions of regular solutions strictly gov- 
erned by the seasonal forcing and more ENSO-like irregular 
solutions, indicating a structural instability of the model. 
In Figure 8 (a) we have plotted such a regular solution ob- 
tained for n = 0.4088, and in Figure 8 (b) an irregular one 
on the other side of the parameter-space discontinuity; for 
Ti = 0.4198. In both of these cases, c = 1, k = 100 and 
d — 1. This irregular solution exhibits a seasonal cycle of 
variable strength interrupted by stronger episodical oscilla- 
tions with some reminiscence of observed ENSO episodes. 
However, by inspection of the time series the similarity to 
the Nino 3 signal is not too convincing. As can be seen from 
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Figure 10. (a): Histogram of signal in Figure 6(a). 
(b): Histogram of signal in Figure 6(b). Conditionally 
averaged evolution of ENSO signal after onset of ENSO 
episode. For simulation with stochastic forcing shown in 
Figure 6(a) (red triangles). For Nino 3 index shown in 
Figure 6(b) (blue squares). The error bars denote the 
standard deviation of the (s n (St)). 
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Figure 6, the solution with stochastic forcing added exhibits 
more qualitative similarity to the Nino 3 signal than the sig- 
nal in Figure 8(b). This qualitative discrepancy is also re- 
flected by the test of determinism. Figure 9(a) shows A(r) 
for the signal in Figure 8(b) and for its phase-randomized 
version. The determinism is higher than the one for the Nino 
3 signal shown in Figure 3(a), and it is less reduced by phase 
randomization. The former is obviously due to the higher 
degree of randomness apparent in the Nino 3 signal. The 
lack of reduction of determinism after phase-randomization 
of the Figure 8(b)-signal is due to the strong seasonal os- 
cillation. Since a major fraction of the power resides in 
the annual Fourier component the phase randomization will 
change the phase of this component, but the signal will still 
have a deterministic appearance after this change. This is 
all quite trivial, so the crucial test is to compute determin- 
ism after de-seasonalization. For the signal in Figure 8(b) 
subtraction of the climatology does not offer an effective fil- 
tering of the seasonal component, since its amplitude varies 
a lot. For the same reason Fourier filtering is also inef- 
fective. In Figure 8(c) we show the de-seasonalized signal 
after application of a mexican-hat wavelet filter. This sig- 
nal represents the "true" ENSO episodes according to the 
time-delay model without stochastic forcing. The computed 
determinism A(r) of this signal is shown in Figure 9(b). The 
determinism is very high, and slows clearly that the ENSO 
episodes in this model are the result of deterministic, low- 
dimcnsional dynamics. The strong reduction of determinism 
after phase-randomization demonstrates that this dynamics 
is nonlinear, and (since the timing of the episodes seems ran- 
dom) chaotic. This is in strong contrast to the results shown 
in Figure 5 for de-seasonalized Nino 3 and signals generated 
by the time-delay model with stochastic forcing, which show 
signals dominated by a stochastic component. 

In the following, we go back to the time-delay simulation 
with stochastic forcing and demonstrate in two more exam- 
ples its similarity to the Nino 3 data. In Figure 10(a) we plot 
a histogram for the Nino 3 data and for the simulation after 
the climatology has been subtracted (the signals in Figure 
6). The ENSO events contribute to the tails of these dis- 
tributions, which are somewhat heavier on the positive side 
due to the relative strength of El Nino compared to La Nina. 
The average time-evolution of ENSO events in Nino 3 and 
simulation can be investigated by means of a superposed- 
epoch analysis (conditional averaging). This analysis works 
as follows: Consider a signal s(t). We define the onsets of 
ENSO events as the times t„, n = 1, . . . , N, for which the 
signal ascends through two standard deviations from the 
mean. Then we produce an ensemble of conditional sig- 
nals Sn(St) — s(t n +5t) and produce the conditional average 
(s„(6t)) = (1/N) s «( St )- Figure 10(c) displays the av- 

erage ENSO structure for Nino 3 (squares) and simulation 
(triangles) computed this way, and shows that the simula- 
tion predicts fairly well the average evolution of ENSO up 
to at least 15 months after the onset of the episode. 

6. Conclusion 

From of a model for ENSO activity similar to the one 
studied here, Tziperman et al. [1995] came to the conclu- 
sion that depending on the strength of the coupling be- 
tween the ocean and the atmosphere, the dynamics could 
undergo quasi-periodicity routes to chaos. The same au- 
thors suggested that the ENSO might be described as a low- 
dimensional chaotic dynamics driven by the seasonal cycle, 
where the appearance of the chaos is due to the nonlinear 
resonance between the natural oscillator of the atmosphere- 
ocean coupling and the seasonal cycle. Further, the irregu- 
lar jumps of the state between different resonances should 



be the cause of the chaos. Unfortunately, SST time series in 
the eastern Pacific are too short for chaos to be proven by 
standard methods. 

In this paper we have applied a simple test for determin- 
ism to show that Nino 3 data, which is the proxy of SST 
in the eastern Pacific, is most adequately described as a 
stochastic process after the seasonal cycle has been removed, 
implying that ENSO is a consequence of mainly stochastic 
and not low-dimensional, chaotic dynamics. A similar con- 
clusion was made by Binder and Wilches [2002], where a 
test for determinism was applied to the Southern Oscilla- 
tion index (SOI) series. SOI measures a pressure difference 
between Tahiti and Darwin in the Pacific Ocean and can 
also be used as a proxy for the El-Nino dynamics. 

Despite the stochastic nature of the Nino 3 signal, we 
have demonstrated here that low-dimensional, deterministic 
dynamics may also be involved. The equatorial- wave equa- 
tion from Tziperman et al. [1994] exhibits such dynamics 
without stochastic forcing, but by adding stochastic forcing 
to this equation and comparing determinism and average 
ENSO structure, with those of Nino 3 data, we conclude 
that this statistical-dynamical model can reproduce impor- 
tant aspects of ENSO dynamics. 
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